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j> \ 1 INTRODUCTION 
•i—i . 

^ , While examining an NRAO VLA Sky Survey (NVSS) image (Con- 
5^ ' don et al. 1998) of the region surrounding the Galactic supernova 
o3 1 remnant G63-7+1-1, an arc of extended emission was noticed near 
RA= 19 ;i 50'"24 s , Dec= 28°16'25" (see Fig. 1). As it resembles 
part of a ring and lies in the Galactic plane, it was identified as 
a possible supernova remnant (SNR). Tian and Leahy 2006 also 
noticed this structure in Canadian Galactic Plane Survey (CGPS; 
Taylor et al. 2003) data with similar sensitivity to but poorer reso- 
lution than the NVSS data, and proposed it as a possible SNR. It 
has now been observed with the VLA and the newly commissioned 
Arcminute Microkelvin Imager (AMI; see AMI Consortium: Zwart 
et al. 2008). These observations are described in Section 2, and cal- 
ibration and data reduction in Section 3. Results are presented in 
Section 4, and we analyse and discuss all available data in Section 
5, confirming G64-5+0-9 as a new Galactic SNR. 



ABSTRACT 

We present observations between 1.4 and 18 GHz confirming that G64-5+0-9 is new Galactic 
shell supernova remnant, using the Very Large Array and the Arcminute Microkelvin Imager. 
The remnant is a shell ~ 8' in diameter with a spectral index of a = 0.47 ± 0.03 (with a 
defined such that flux density S varies with frequency vasS« v _a ). There is also emission 
near the centre of the shell, ~ 1' in extent, with a spectral index of a = 0.81 ±0.02. We do not 
find any evidence for spectral breaks for either source within our frequency range. The nature 
of the central object is unclear and requires further investigation, but we argue that is most 
unlikely to be extragalactic. It is difficult to avoid the conclusion that it is associated with the 
shell, although its spectrum is very unlike that of known pulsar wind nebulae. 
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Figure 1. NVSS image of G64-5+0-9; the resolution is 45 arcseconds and 
the map noise is 3 mjy beam - 1 . 



* We request that any reference to this paper cites "AMI Consortium: 
Hurley- Walker et al. 2009" 

t Issuing author; E-mail: nh313@mrao.cam.ac.uk 



2 OBSERVATIONS 

The position of the source close on the sky to the ring centre was 
chosen as the pointing centre for observations with the VLA and 
AMI. AMI is a dual set of interferometric arrays located at the 
Mullard Radio Astronomy Observatory, Lord's Bridge, Cambridge, 
UK. The AMI Small Array (SA) consists of ten 3.7-m diameter 
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Table 1. Observations of G64-5+0-9. 



Telescope 


AMI SA 


AMI LA 




VLA 




Date 


2007 October 


2008 October 




2008 July 




Flux Calibrator 


3C286 


3C48 




3C48 




Phase Calibrator 


J2023+3153 


J2023+3153 




J 1925+2 11 




Frequency/ GHz 


14-18 


14-18 


1.43 


4.86 


8.46 


Observation(s) Length 


16 hours 


10 hours 


40 minutes 


40 minutes 


20 minutes 


Synthesised beam FWHM 


3-2' 


1-0.3' 


0.75' 


0.22' 


0.12' 


Primary beam FWHM 


22-16' 


5.6-5.0' 


32' 


9.3' 


5.3' 



equatorially-mounted dishes with a baseline range of ~ 5 - 20 m, 
while the AMI Large Array (LA) has eight 12.8-m-diameter dishes 
with a baseline range of ~ 20 - 100 m. Both arrays observe I+Q in 
the band 12-18 GHz, each with system temperatures of about 25 K. 

The backends are analogue Fourier transform spectrometers, 
from which the complex signals in each of eight channels of 750- 
MHz bandwidth are synthesised, and the signals in the synthesised 
channels are correlated at the ~ 10 per cent level. In practice, the 
lowest two frequency channels are generally not used due to a poor 
correlator response in this frequency range and, currently, interfer- 
ence. 

Frequencies, resolutions, dates, calibrators and observing 
times for each observation are shown in Table 1. The AMI SA 
observed a single pointing while the AMI LA used a seven-point 
raster mode with 4' spacing. J2023+3153 (RA= 20 /l 23 m 19 s .O, 
Dec= 31°53'02".3) was chosen as a phase calibrator from the Jo- 
drell Bank VLA Survey (JVAS; Patnaik et al. 1992; Browne et al. 
1998; Wilkinson et al. 1998) on the basis of its proximity and flux 
density (1.1 Jy at 15 GHz). 

During the time of our VLA observations, about half of the an- 
tennas of the VLA were upgraded Extended VLA (EVLA) anten- 
nas. This necessitated an extra step of calibration on VLA-EVLA 
baselines. The telescope was in D-array, and we used the L, C and 
X band receivers at centre frequencies of 1.43, 4.86 and 8.46 GHz 
respectively. J1925+211 (RA= 19 ,! 25 m 59' 5 .6, Dec= 21°06'26".2) 
was selected from the VLA Calibrator Manual 1 as the phase cali- 
brator at 1.43, 4.86 and 8.46 GHz in D-array; it has an assumed flux 
density of 1.3, 1.5 and 1.0 Jy in each of these bands, respectively. 



3 CALIBRATION AND DATA REDUCTION 

The AMI data reduction was performed using our software tool 
REDUCE. This is used to apply path-compensator and path-delay 
corrections, to flag interference, shadowing and hardware errors, to 
apply phase and amplitude calibrations and to Fourier transform the 
correlator data readout to synthesise the frequency channels, before 
outputting to disk in wv-FITS format suitable for imaging in AIPS. 

Flux calibration was performed using short observations of 
3C48 and 3C286 near the beginning and end of each run, with as- 
sumed I + Q flux densities for these sources in the AMI channels 
consistent with Baars et al. (1977) (see Table 2). As Baars et al. 
measure I and AMI measures I + Q, these flux densities include 
corrections for the polarization of the sources derived by interpo- 
lating from VLA 5-, 8- and 22-GHz observations. After phase cal- 
ibration, the phase of AMI over one hour is generally stable to 5° 
for channels 4-7, and to 10° for channels 3 and 8. 

The system temperatures of each AMI antenna are continously 



Table 2. Assumed I + Q flux densities of 3C286 and 3C48 over the 
commonly-used AMI band. 



Channel 


v/GHz 


S icm /Jy 




3 


14.2 


3.61 


1.73 


4 


15.0 


3.49 


1.65 


5 


15.7 


3.37 


1.57 


6 


16.4 


3.26 


1.49 


7 


17.1 


3.16 


1.43 


8 


17.9 


3.06 


1.37 



monitored using a modulated noise signal injected at each an- 
tenna; this is used to continuously correct the amplitude scale in 
a frequency-independent way. The overall consistency of the flux- 
density scale is estimated to be better than 5 per cent. 

Since the AMI antennas are sensitive to I+Q and are equa- 
torially mounted, this polarization is fixed on the sky during the 
observation. Q may be positive or negative, but is expected to be 
small when integrated over the whole source. 

VLA data reduction was performed entirely within AIPS. The 
VLA-EVLA baselines were also calibrated using BLCAL as de- 
scribed in the guidelines for post-processing EVLA data in AIPS 2 . 
This corrects for closure errors on these baselines caused by non- 
matched EVLA and VLA bandpass shapes. 

Maps were made using IMAGR in AIPS from each channel of 
the AMI SA, AMI LA and the VLA (Fig. 2). Combined-channel 
maps of the AMI SA and AMI LA observations are shown in Fig. 
3 and Fig. 4, respectively. The VLA and AMI SA images shown 
are not corrected for the primary beams of the telescopes, but when 
measuring flux densities this correction is made. In order to pro- 
duce a raster map, the AMI LA images are of course primary-beam 
corrected. The signal-to-noise of the 1.43-GHz image was large 
enough that self-calibration from the first non-negative CLEAN 
components was used to enhance the image. 



4 RESULTS 

We note that the maps show two sources of emission: a ring-like 
structure surrounding a slightly extended, apparently central object. 
In those maps that do not resolve the ring completely, the emission 
from these two sources may be confused. 

As we are using interferometers sensitive to a limited range 
of angular scales, it is not possible to obtain accurate flux densities 
for structures on all scales from all of the observations. Only the 
VLA 1.43-GHz and AMI SA observations have the required short 
baselines to provide sensitivity to the extended ring emission. The 
AMI SA does not have the required resolution to measure the flux 
density of the central emission without confusion with the ring. 

To measure flux density, we adopt the fitting method of Green 



http://www.aoc.nrao.edu/~gtaylor/csource.html 



2 http://www.vla.nrao.edu/astro/guides/evlareturn/postproc.shtml 
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Figure 2. VLA images of G64-5+0-9. Top: 1.43 GHz, a rms = 
110 jUjy beam -1 , uniform weighting; middle: 4.86 GHz, a rms = 
38.5 lljy beam -1 , uniform weighting; bottom: 8.46GHz, a rms = 
43.5 fiJy beam -1 , natural weighting with three-sigma contours. The box 
in the middle panel indicates the area that is shown in the lower panel. In 
this and subsequent images, the FWHM of the CLEAN restoring beam of 
the observation is shown as an ellipse inside a box on the lower left. 



(2007), using the program FITFLUX. In this method, a flux den- 
sity is fitted by drawing a polygon around the object and fitting a 
tilted plane to the pixels around the edges of the polygon. The tilted 
plane is then removed from the image before integrating the emis- 
sion within the polygon. This method is sensitive to the fitting area 
selected but only at the one per cent level, which is included in the 
estimate of the error. This method allows us to remove contribution 
from the background, especially the contribution from the ring in 
the case of measuring the flux density of the central emission. 

The main sources of error on the flux density measurements 
are the thermal noise Cmis and the error on the flux calibration 
Os- Long-term measurements of the AMI flux calibration show it 
to have an r.m.s. error of three per cent. The VLA flux calibra- 
tion is performed using GETJY which scales the fluxes of the field 
and phase calibrator using the flux density of 3C48 consistent with 
Baars et al.. This process has an error of less than one per cent for 
each of the frequency channels. Therefore to obtain the error on a 
flux density measurement, the thermal noise of each map is added 
in quadrature with the flux calibration error of the telescope used, 
and a one per cent error when FITFLUX is used. 

We now discuss the results from each observation in turn. 



4.1 VLA 1.43 GHz 

Fig. 2 shows that this object is clearly a shell SNR, brighter in the 
north-east. Interestingly, the central emission is also brightest in 
this direction. The full ring structure is clearly discernable in this 
map, as compared to the partial structure shown in the NVSS map 
(Fig. 1). This is due to the enhanced uv coverage and lower noise 
of the longer pointed observation. 

FITFLUX was used to obtain the flux density of the central ob- 
ject as 17. 1 ± 0.8 mjy, and the flux density of the entire remnant as 
1 19 ± 5 mjy. The former was subtracted from the latter to obtain 
the flux density of the ring as 102 ± 5 mjy. 



4.2 VLA 4.85 GHz 

At 4.85 GHz, the ring is resolved out to the extent that it is not 
possible to extract a reliable estimate of its integrated flux density. 
However, the structure of the emission close to the ring centre be- 
comes more clear; its integrated flux density is listed in Table 3. 

The north-east section of the ring shows possible contamina- 
tion with a point source. Using the MPS routine SLICE, a cross- 
sectional profile across this source, perpendicular to a tangent along 
the ring at this point, was produced in order to measure the flux 
levels. The base flux level from the ring is 220 /iJy beam -1 and 
the peak point source flux is 620 /xJy beam -1 . The width of the 
source profile at 420 /iJy beam -1 is 12", which is the resolution of 
the telescope at this frequency. Therefore the source is unresolved 
and we postulate it is a background contaminating source. FITFLUX 
gives a flux density estimate of 0.41 ± 0.04 mjy for this source. 



4.3 VLA 8.46 GHz 

At 8.46 GHz, the ring is almost entirely resolved out. However we 
obtain good information on the structure of the near-central emis- 
sion, as shown in the lowest panel of Fig. 2. FITFLUX gives a mea- 
surement of the flux density of this source, shown in Table 3. 
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Figure 3. AMI SA combined-channel (14-18 GHz) image of G64-5+0-9; 
the map noise is 192 /iJy beam -1 . 




19 h 51 m 00 s 50 m 45 s 50 m 30 s 50 m 15 s 50™00 s 



Figure 4. Map of a primary-beam-corrected combined-channel (14- 
18 GHz) seven-point raster of G64-5+0-9 observed by the AMI LA. The 
central pointing has a thermal noise of 32 /iJy beam -1 and the surrounding 
six pointings have a thermal noise of 58 /iJy beam -1 . 



4.4 AMI 14-18 GHz 

The AMI SA, due to its sensitivity to larger angular scales, picks 
up much more large-scale Galactic emission in the region (see Sec- 
tion 4.5 and Fig. 3). The north-east part of the ring structure is clear 
but is confused with the central source. 

The AMI LA combined-channel raster map is shown in Fig. 4. 
The ring is slightly resolved out by the AMI LA, but the source ap- 
parently near the centre shows up clearly and has similar structure 
to that in the VLA maps at lower frequency. FITFLUX was used to 
find the flux density of this near-central emission for each channel 
and these results are listed in Table 3. It is possible that the AMI 
LA is slightly resolving out this object at the higher end of its fre- 
quency coverage. However the slight drop in flux density is within 
the error bars. 

FITFLUX was used to find the flux density of the whole object 
from the AMI SA combined-channel map, and the flux density of 
the central source from the AMI LA combined-channel map. These 
were 35.3 ± 1.4 mjy and 2.49 ± O.lOmJy, respectively. The latter 



Table 3. VLA and AMI LA I+Q flux densities of the emission close to the 



ring centre. 


v/GHz 


Sj/mJy 


1.43 


17.07 ±0.36 


4.86 


7.92 ±0.16 


8.46 


4.01 ±0.08 


14.2 


2.69 ±0.17 


15.0 


2.58 ±0.12 


15.7 


2.50 ±0.1 3 


16.4 


2.48 ±0.12 


17.1 


2.27 ±0.12 


17.9 


2.02 ±0.11 



was subtracted from the former to produce an estimate of the flux 
density of the ring at 16 GHz of 32.8 ± 1 .4 mjy. 

4.5 Data from the literature 

Most Galactic radio surveys do not possess the required resolution 
to resolve G64-5+0-9 well, as it is only 8' in diameter. The Effels- 
berg 2.7-GHz survey data (Reich et al. 1984) show a small knot of 
emission at the SNR's location, and a large amount of emission on 
broader angular scales. It was not possible to extract meaningful 
flux densities for the ring and central source from these data. 

This object is also covered by the Westerbork Synthesis Ra- 
dio Telescope Galactic plane survey 3 (see e.g. Taylor et al. 1996) 
at 327 MHz, which provides low-frequency, large-scale informa- 
tion, at a resolution of l'x2'.ll at this declination. Unfortunately 
the residual grating ring of a bright, distant source lies directly over 
the objects of interest. 

Infrared data from the Infrared Astronomical Satellite 
(IRAS;NASA RP-1190 1987) are too low in resolution to provide 
any useful information about the region. The Second Digitised Sky 
Survey (DSS2) 4 shows only foreground stars in the Galaxy. The 
ATNF Pulsar Catalogue 5 (see Manchester 2005) shows no pulsars 
within one degree of the remnant. 



5 DISCUSSION 

For both the ring and the near-central emission we fit power-law 
spectra. We define a such that flux density S varies with frequency 
V as5« v - ". Spectra were fitted to these data using a Gaussian 
likelihood function sampled by a Markov Chain Monte Carlo tech- 
nique. This method copes with asymmetric errors in log 5 across the 
frequency channels and provides an error estimate on the spectral 
index directly from the posterior distribution. 

Using the VLA 1 .43 GHz and AMI measurements of the 
flux density of the ring, we calculate that it has a spectrum with 
a = 0.47 ±0.03. The structure and spectral index identify the shell 
source G64-5+0-9 as a shell SNR. The surface brightness of the 
SNR at 1 GHz is ~ 3 x 10 -22 W m -2 Hz -1 sr -1 , which is faint for 
known SNRs (in the faintest 10 per cent of identified SNRs - see 
e.g. Green 2004). This relatively low surface brightness suggests 
that the remnant is old, but on the other hand, its highly circular 
shell suggests it is young. (Tian and Leahy 2006) suggest a possible 
distance of 1 1 kpc, based on possibly related HI features, at which 

3 http://www.ras.ucalgary.ca/wsrt_survey.html 

4 http://archive.eso.org/dss/dss 

5 http://www.atnf.csiro.au/research/pulsar/psrcat/ 
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it would have a diameter of ~26 pc. This is not the first SNR visi- 
ble in NVSS data which has been overlooked (see e.g. G353-9— 2-0 
identified by Green 2001). 

Fitting a spectrum to the data in Table 3, we find that the cen- 
tral emission has a steep spectrum with a = 0.81 ±0.02 (Fig. 5). 
This emission is of interest, because of both its apparent position 
relative to the ring and its structure. It could be a background, un- 
related, extragalactic source, or a Galactic source that may or may 
not be related to the shell remnant. 

First we consider a source at redshift z around 0.1. Then its 
angular size of around 1' implies a physical size D of 100 kpc and 
its luminosity P\ .43 is around 3 x 10 22 WHz _1 sr _1 (taking Ho = 
72kms~ 1 Mpc~ 1 ). These values are typical of FR-I (Fanaroff and 
Riley 1974) twin-jet or tail radio sources powered by jet-producing 
machines. However the emission is not reminiscent of known twin- 
jets or tails - particularly at 8 GHz - and the lack of spectral steep- 
ening from 1.43 to 18 GHz also tells against it being this type of ra- 
dio source. A substantially higher redshift for the emission appears 
ruled out: its structure is very difficult to reconcile with that of a 
powerful radio source of whatever physical size or orientation, and 
the emission shows no sign of the synchrotron or inverse-Compton 
losses associated with the extended regions of powerful sources. 

'Halo' emission from certain galaxy clusters might resem- 
ble the shape of the near-central emission. However, 'halo' emis- 
sion has an extremely steep radio spectrum (see e.g. Brunetti et al. 
2008). 

The remaining extragalactic possibility would seem to be syn- 
chrotron emission from a nearby system of interacting galaxies. At 
a redshift z —0.02, its angular size of around 1' implies a physical 
size D of 20 kpc and its flux density Si 43GHz of 17.1 mjy implies 
a luminosity .Pi.43GHz of 1 x 10 21 WHz _1 sr _1 Thus, this source's 
angular size, flux density and structure are consistent with syn- 
chrotron emission from two or three interacting galaxies at z of 
around 0.01 - 0.05. 

However, from NVSS source counts 6 , we estimate at 15 per 
cent the probability of a source of Si 43 ^17mJy appearing to lie 
inside the shell of the SNR purely by chance. Thus the probability 
of any extragalactic radio source of this flux density and of such 
large angular extent, lying within one arcminute of the centre, is 
minimal. 

Turning to possible Galactic sources of emission for the near- 
central source, chance assocation along the line-of-sight with ther- 
mal Galactic sources such as HII regions is excluded both by the 
non-thermal radio spectrum of the central source, and its lack of 
infra-red emission. The location of the extended emission near the 
centre of the shell suggests that G64-5+0-9 may be a 'composite' 
SNR, i.e. a shell containing a central pulsar wind nebula (PWN) 
- for example, see G328-1-1-8 (Dickel et al. 2000). However, 
the radio spectral index of ~ 0.8 for the central source is con- 
siderably larger than expected, since pulsar wind nebulae usually 
have relatively flat spectral indices at gigahertz frequencies (e.g. 
Kothes 2008). If the central source is a pulsar wind nebula asso- 
ciated with G64-5+0-9, then it requires a spectral break at below 
1.4 GHz, which would be unusual. A low frequency spectral break, 
at about 1 .3 GHz - three times lower than for any other PWN - has 
been identified in DA 495 (=G65-7+l-2, Kothes et al.). DA 495 
is thought to be rather old (~ 20, 000 years), whereas the highly 
circular nature of the shell of G64-5+0-9 suggests it is relatively 
young. 



6 http://sundog.stsci.edu/first/catalog_paper/node6.html 
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Figure 5. Logarithmic plot of the flux density measurements of the source 
seen near the centre of G64-5+0-9 in 3. The dashed line indicates a fitted 
spectrum of a = 0.81 using the method described in Section 5. 

To identify the nature of this object and fully constrain its as- 
sociation with the remnant, HI absorption-line observations are vi- 
tal. 



6 CONCLUSIONS 

Follow-up obervations of an arc of emission seen in NVSS and 
CGPS data have been made at frequencies from 1.4 to 18 GHz. We 
draw the following conclusions from our observations: 

• G64-5+0-9 is confirmed as a shell SNR and its emission fol- 
lows a power-law spectrum with a = 0.47 ± 0.03; 

• by considering its size, shape and spectral behaviour, and a 
probabilistic analysis using NVSS radio source counts, we argue 
that the emission apparently at the centre of the ring is most un- 
likely to be any kind of extragalactic radio source; 

• it is difficult to avoid the conclusion that the apparently cen- 
tral emission is associated with the shell, although its spectrum is 
very unlike that of a pulsar wind nebula, and it therefore merits 
further investigation, particularly follow-up HI absorption line ob- 
servations. 
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